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Turbulent flame propagation 
in partially premixed flames 

By T. Poinsot 1 , D. Veynante 2 , A. Trouve 3 4 AND G. Ruetsch 3 


1. Introduction 

Turbulent premixed flame propagation is essential in many practical devices. In 
the past, fundamental and modeling studies of propagating flames have generally 
focused on turbulent flame propagation in mixtures of homogeneous composition, 
i.e. a mixture where the fuel-oxidizer mass ratio, or equivalence ratio, is uniform. 
This situation corresponds to the ideal case of perfect premixing between fuel and 
oxidizer. In practical situations, however, deviations from this ideal case occur 
frequently. In stratified reciprocating engines, fuel injection and large-scale flow 
motions are fine-tuned to create a mean gradient of equivalence ratio in the com- 
bustion chamber which provides additional control on combustion performance. In 
aircraft engines, combustion occurs with fuel and secondary air injected at, vari- 
ous locations resulting in a nonuniform equivalence ratio. In both examples, mean 
values of the equivalence ratio can exhibit strong spatial and temporal variations. 
These variations in mixture composition are particularly significant in engines that 
use direct fuel injection into the combustion chamber. In this case, the liquid fuel 
does not always completely vaporize and mix before combustion occurs, resulting 
in persistent rich and lean pockets into which the turbulent flame propagates. 

From a practical point of view, there are several basic and important issues re- 
garding partially premixed combustion that need to be resolved. Two such issues 
are how reactant composition inhomogeneities affect the laminar and turbulent 
flame speeds, and how the burnt gas temperature varies as a function of these in- 
homogeneities. Knowledge of the flame speed is critical in optimizing combustion 
performance, and the minimization of pollutant emissions relies heavily on the tem- 
perature in the burnt gases. Another application of partially premixed combustion is 
found in the field of active control of turbulent combustion. One possible technique 
of active control consists of pulsating the fuel flow rate and thereby modulating the 
equivalence ratio (Bloxsidge et al 1987). Models of partially premixed combustion 
would be extremely useful in addressing all these questions related to practical sys- 
tems. Unfortunately, the lack of a fundamental understanding regarding partially 
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premixed combustion has resulted in an absence of models which accurately capture 
the complex nature of these flames. 

Previous work on partially premixed combustion has focused primarily on laminar 
triple flames. Triple flames correspond to an extreme case where fuel and oxidizer 
are initially totally separated (Veynante et al. 1994 and Ruetsch et a l 1995). 
These flames have a nontrivial propagation speed and are believed to be a key 
element in the stabilization process of jet diffusion flames. Different theories have 
also been proposed in the literature to describe a turbulent flame propagating in 
a mixture with variable equivalence ratio (Muller et al. 1994), but few validations 
are available. The objective of the present study is to provide basic information 
on the effects of partial premixing in turbulent combustion. In the following, we 
use direct numerical simulations to study laminar and turbulent flame propagation 
with variable equivalence ratio. 

2. Framework for analyzing and modeling partially premixed flames 

Perfectly premixed combustion is usually described using a progress variable 
c = (Y$ — Yf)/(Yf — Yp), where Yp and Yo are the fuel and oxidizer mass frac- 
tions, and the superscripts 0 and 1 refer to the values in the unburnt reactants and 
burnt products, respectively. Using the assumption of single-step chemistry and 
unity Lewis numbers, the progress variable provides a complete description of the 
transition from unburnt to burnt states and is the single relevant quantity used 
in model development and the postprocessing of simulation results. In partially 
premixed combustion, a new theoretical framework is required which will allow 
variable equivalence ratio along with simultaneous premixed and diffusion modes of 
combustion. This framework must use at least two scalar variables: one variable to 
describe the species composition, and a second variable to describe the progress of 
the premixed reaction. We use the mixture fraction Z as a description of the species 
composition, and a modified form of the progress variable c which accommodates 
the variable species composition in the fresh reactants. 

We assume irreversible single step chemistry and unity Lewis numbers: 

F + r s (0 + bN 2 )^ P (1) 

where r 9 is the stoichiometric oxidizer-fuel mass ratio and b is the N 2 -0 mass ratio, 
where N 2 is a diluent in the fresh reactants. The mixture fraction Z is then defined 
as: 

Z = ~ Y o/ r » + VIM 1 + *)) (0) 

l + l/(r 9 (l + b)) * 

For stoichiometric mixtures, Z is equal to Z 9i — l/(r 9 (l + b) + 1).* The fuel 

* The equivalence ratio, (f) = V 9 \ p jY(), is a more familiar quantity to the engineering com- 
munity. However, <j) and Z are simply related (Muller et al. 1994) through <f> = Z(1 — 
Z 9 t)/Z s t/( 1 — Z). We use Z rather than <p for the following reasons: (f> is a conditional quan- 
tity that is only defined in the unburnt reactants whereas Z is not only defined everywhere in the 
flow, but is also conserved through the reaction, and Z is a linear combination of the species mass 
fractions and leads to straightforward expressions when averaging is employed. 
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FIGURE 1. Fuel and oxidizer mass fractions as a function of the mixture fraction 

z. : Fuel Yp/Yj? (mixing line) ; : Oxidizer Yg/Y'ftft (mixing line) ; 

: Fuel Y'p/Y'f 1 (fast chemistry line) ; : Oxidizer Yq/Y'q° (fast chemistry 

line). The arrows indicate the transition from unburnt to burnt states in the case 
of perfectly premixed combustion, at a given equivalence ratio. 


consumption rate, wp, and the heat release rate, wp, are written as: 

Cjf = Ap v YpY'o exp(-T a /T) and thr = q^F (3) 


where .4, p , n, m are model constants, p is the mass density, and q is the heat 
of reaction per unit mass of fuel. The activation temperature T a is specified via 
a Zeldovich number, ft = aT a /T b (Z st ), where a is the stoichiometric heat release 
factor, q = (T b (Z 3t ) ~ T 0 )/T b {Z 9t ). The unburnt gas temperature, T 0 , is assumed 
uniform in the present study, and the adiabatic flame temperature, Tb{Z st ), is 
calculated under stoichiometric conditions. 

The mixture composition upstream of the flame zone is only a function of Z and 
may be described by the two mixing lines shown in Fig. 1 and given by 


yO 

— Z and 

t'OO 

1 F 


y ° 

l°- = l-z 

V3C 

x o 


(4) 


where the superscript oo denotes the value of the mass fraction in the respective 
feeding streams, such that Yp? = 1 and Yq 1 = 1/(1 + b). 

If the chemistry is sufficiently fast, the mixture composition downstream of the 
premixed flame corresponds to the classical Burke-Schumann limit where fuel and 
oxidizer cannot coexist. This limit is a function of Z alone, and is sketched in Fig. 1 
and given by 

^ -*« («•££) “ d ,5) 

Premixed combustion changes the mixture composition from an unburnt state, as 
described by Eq. (4), to a burnt state, as described by Eq. (5). Under the flamelet 
assumption, this change occurs in a thin flame zone. Note that in perfectly premixed 
combustion, the mixture fraction Z is constant and the transition from unburnt to 
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burnt states occurs on a vertical line in Fig. 1. In partially premixed combustion, 
the transition may occur with simultaneous variations in Z. 

Equations (4) and (5) lead to the following generalized definition of the premixed 
reaction progress variable: 


ZYf - Y f 


ZYp° - Max 


(o. tt) 

For lean mixtures, where Z < Z st everywhere, we have 

Y f 


vx 
1 F 


C — 1 — 


ZY? 


and for rich mixtures, with Z > Z s , everywhere, the following holds: 

Yo 


ZYf - Yf 


ZYj? 


-(f 5fc) 


= 1 - 


yx 
1 F 


(1 -Z)Yg> 


(6) 

(7) 

( 8 ) 


If Z is constant, Eq. (6) reduces to the standard definition of c used in perfectly 
premixed combustion. 

For the sake of simplicity, we now limit our discussion to the case of a lean mixture. 
A balance equation for c may be derived from basic conservation equations for the 
fuel mass fraction Yf and for the mixture fraction Z\ 


dc dc _ \ d t dc\ Uf 2 D dc dZ 

dt +U, !h r i ~ ~pdT k \ pD d^ k ) ~ pZYp + ~Yd7idxi (9) 

where u l is the fluid velocity and D is the mass diffusivity. This equation is similar to 
the one obtained in perfectly premixed combustion, except for the last term on the 
right-hand side. The sign of this additional term can be either positive or negative, 
suggesting flame propagation can either accelerate or decelerate as a result of partial 
premixing. Following Trouve and Poinsot (1994), the conservation equation for c 
may be used to define the displacement speed of iso-c surface contours: 


w(c — c * ) = 


|Vc| 


dc 


dc 


dt +u ‘dx i 


|Vc| 


~V • (pDVc) 

P pZYZ? 


2D 

+ ~Z Vc ' 


vz 


where all quantities are evaluated at c = c* . 
is: 


(10) 

An alternative form of this equation 


w 


I Vc j [p 


-V • (pDVc) 


u>F 

pzyf 


2D 


n VZ 


(ID 


where n is the local unit vector normal to the iso-c surface, n = — Vc/ |Vc|. Adopt- 
ing a flamelet point of view, we identify the thin flame surface as an iso-c surface 
with c* = 0.8. Equation (11) can then be interpreted as an expression for the flame 
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propagation speed. The terms within brackets on the right-hand side of Eq. (11) 
show the dependence of the flame propagation speed on the local mixture frac- 
tion Z. The last term on the right-hand side shows the dependence of the flame 
propagation speed on the local Z-gradient normal to the flame. Hence, one basic 
effect of incomplete reactant mixing is the modification of the local flame speed, 
w(Z, n • VZ). 

We now discuss the implications of partial premixing in the framework of flamelet 
combustion. In the flamelet picture, the mean reaction rate may be written as the 
product of a mean mass burning rate times the flame surface density: 

(wf) = <m) s E (12) 

where in is the local mass burning rate per unit flame surface area, in = f n d >rdn, 
and E is the mean flame surface-to- volume ratio (the flame surface density). The 
operator () 5 denotes a flame surface average (Pope 1988). 

Partial premixing can induce modifications of the mean reaction rate through 
several mechanisms: a modification of the local flame structure and corresponding 
modifications to the mean mass burning rate (m) s , and contributions to the flame 
wrinkling resulting in a modification to the flame surface density E. The effect 
of partial premixing on flame wrinkling may be analyzed by considering the exact 
balance equation for E (Pope 1988, Candel &; Poinsot 1990, Trouve & Poinsot 1994): 


^ + V-(u) s E + V-(u’n) s E = (k) s E = (V * u - nn : S7u) s E + (u>V • n) s E (13) 

where k is the flame stretch, which is decomposed in Eq. (13) into a production 
term due to hydrodynamic straining and a production or dissipation term due to 
flame propagation. The propagation term is the mean product of the local flame 
propagation speed, te, times the local flame surface curvature, V • n. Hence, the 
effect of partial premixing on the local flame speed, te(Z, n-VZ), as seen in Eq. (11), 
can be interpreted as an effect of partial premixing on flame stretch, k(Z, n • VZ), 
and thereby an effect on E. One objective of the present study is to determine the 
relative weight of effects induced by partial premixing on E and (m) s relative to 
the effects of turbulence on these quantities. 

3. Numerical configurations and diagnostics 

In the present study, one-, two-, and three-dimensional direct numerical simula- 
tions are performed with variable density and simple chemistry. The simulations 
use a modified Pade scheme for spatial differentiation that is sixth-order accurate 
(Lele 1992), a third-order Runge-Kutta method for temporal differentiation, and 
boundary conditions specified with the Navier Stokes characteristic boundary con- 
dition procedure (Poinsot & Lele 1992). We refer the reader to the Proceedings of 
the 1990, 1992, and 1994 CTR Summer Programs for further details concerning the 
system of equations solved and the numerical methods. 
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FIGURE 2. Configurations for DNS of partially premixed flames 


The numerical configuration corresponds to a premixed flame propagating into 
a mixture with variable equivalence ratio. The mixture composition upstream of 
the flame is specified according to the probability density function, an integral 
length scale of the scalar field, and the relevant directions of inhomogeneity. The 
probability density function of Z is denoted as p(Z) and can be characterized by its 
mean and rms values, (Z) and Z'. The amplitude of the fluctuations, A Z, can be 
used in place of the rms fluctuation for laminar cases. The characteristic integral 
length scale of the Z -field is denoted as 1%. The directions of inhomogeneity in the 
Z-field are compared to the mean flame front orientation. We choose the direction 
of mean flame propagation as the x-direction. A numerical configuration is called 
x - i n homogeneou s if gradients of Z exist in the x-direction, the direction normal 
to the mean flame front. Likewise, a y — inhomogeneous configuration corresponds 
to a case where species gradients exist tangent to the mean flame front. A fully 
turbulent three-dimensional configuration is called xyz- inhomogeneous. 

Four different configurations are pursued in this study, as depicted in Fig. 2: Case 
1 is a one- dimensional, x-inhomogeneous, unsteady, laminar flame, with a double- 
peak Z-pdf; Case 2 is a two-dimensional, y- inhomogeneous, steady, laminar flame. 
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with a double-peak Z-pdf; Case 3 is a two-dimensional, zy-inhomogeneous, un- 
steady, laminar flame, with a triple-peak Z-pdf; and Case 4 is a three-dimensional, 
xyz-inhomogeneous, non- stationary, turbulent flame, with a Guassian Z-pdf. These 
configurations each have a slightly different simple chemistry scheme, as summa- 
rized in Table I. Case 4 corresponds to a single step reaction mechanism proposed 
by Westbrook k Dryer (1981) for C 3 H 8 -aiT combustion. 

Table I. Parameters for the four simulation configurations. 


Case 

Dim 

b 

r s 

%st 

P 

n 

rn 

P 

a 

1 

ID 

0 

i 

0.5 

1 

i 

i 

8 

0.75 

2 

2D 

0 

i 

0.5 

2 

i 

i 

8 

0.75 

3 

2D 

0 

i 

0.5 

1 

i 

i 

8 

0.75 

4 

3D 

3.29 

3.64 

0.06 

1.75 

0.1 

1.65 

8 

0.75 


In all cases we characterize the effects of partial premixing by comparing the 
results to those obtained with perfect premixing in the same configuration. The 
effects of partial premixing on the local flame structure are characterized by the 
mass burning rates: 


m 


and 


m 


r = 


(14) 


m{Z st ) m((Z)) 

where the unprimed quantity uses the stoichiometric homogeneous laminar flame 
as a reference, whereas the primed quantity uses the homogeneous laminar flame 
with Z = (Z). The global effects of partial premixing are characterized by the total 


reaction rate ratios: 


R = 


Q 

Do(Z»t) 


and R' = 


Q 


M(Z)) 


(15) 


where 0. = f y (uijr)dV, with Do corresponding to a homogeneous, planar, laminar 
flame. These total reaction rate ratios can be rewritten as: 


WsZdV (Sy) 
f Ed V So 


and 


„ /(r > )sSdV(5y) 

fEdV S 0 


(16) 


where Sv is the total flame surface area within V , and 5 0 is the projected area of 
the flame on a surface perpendicular to the direction of mean propagation. The 
first ratio in these expressions for R and R! accounts for modifications of the mean 
mass burning rate due to partial premixing, and the second ratio accounts for 
flame surface wrinkling due to turbulence and partial premixing. We write W = 
Sv/So and P = f (r') s EdV / J EdV . The effects of partial premixing on flame 
temperatures are characterized by the following temperature ratios: 


6 = 


Tmaz ~ To 

T m az,o(Z 3 l) ~ To 


and 8 ' = 


Tmax — To 
T m ax,o{{Z)) ~ Tq 


(17) 


where T ma i,o corresponds to the case of a homogeneous, planar, laminar flame. 
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FIGURE 3. Case 1: One-dimensional, .r-inhomogeneous, unsteady, laminar flames 


4. Case Is One-dimensional, ar-inhomogeneous, unsteady, laminar flames 

In this case, the inhomogeneity of reactant species is longitudinal with respect 
to the flow, and thus the flame response to temporal fluctuations in the mixture 
composition is studied. The mixture composition is forced at the inlet of the compu- 
tational domain in order to generate harmonic perturbations in the Z-field upstream 
of the flame. These perturbations in Z are characterized by their mean value, (Z), 
their amplitude, A Z, and their wavelength lz (see Fig. 3). This case is well-suited 
to bring basic information on both flame structure modification and quenching by 
partial premixing. 

In order to study the flammability limits of partially premixed flames, it is im- 
portant to determine whether the simplified kinetic scheme used in the simulation 
is capable of reproducing realistic variations of the laminar flame speed, Sl, when 
variations occur in the mixture composition or equivalence ratio, <f>. In particular, 
the lean and rich flammability limits must be correctly predicted. The single-step 
chemistry model presented in Section 2 does not have this capability unless heat 
losses are added to the energy equation. The choice of a nonadiabatic flame may 
be viewed as a simple fix to produce realistic variations S /,(</>), which is presented 
in Fig. 4. Following W illiams (1985), we use a volumetric heat loss term C that is 
linear in (T — To) (see also Poinsot et al. 1991): 


C = 



Dth 


CpTo 


a 

l—o 


(18) 


where h is a model constant, chosen as h = 0.031, and t = (T - T 0 ) / (Tt,(Z st ) - To). 

As seen in Fig. 4, no abrupt transition to extinction is observed for the adiabatic- 
single-step chemistry model, where very lean and very rich mixtures continue to 
burn. As a result, flame speeds are unrealistically high in these lean and rich 
regions. However, a domain of flammability is obtained using nonadiabatic single- 
step chemistry. This domain compares reasonably well to computations performed 
with a detailed mechanism proposed by Westbrook &z Dryer (1981), for CH 4 - air 
flames. While the prediction of the rich flammability limit is overestimated, the 
overall level of accuracy is deemed acceptable at the present, stage. 
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FIGURE 4. Variations of normalized flame speed with equivalence ratio, s(4>), for 
a one-dimensional, homogeneous, laminar flame. □ : Adiabatic one-step chemistry; 
* : Non-adiabatic one-step chemistry; o : Detailed mechanism. 



x (Distance to inlet) 

FIGURE 5. Case 1 with ( Z ) = Zst = 0.5, A Z = 0.2 (at inlet), and lz/^t — 14, 
where 6° L is the thermal thickness of the perfectly premixed flame with Z = (Z). 
+ : Reduced Yf ; : Reduced Yq ; : Reduced temperature r. 


Figure 5 presents a typical snapshot of Vo, and temperature profiles across 
the flame zone. Species mass fractions are normalized in this figure by their stoi- 
chiometric values. One difficulty in these low Reynolds number simulations is that 
the perturbations in Z imposed at the inlet are strongly affected by molecular dif- 
fusion and are significantly damped before they reach the flame. In this situation, 
the flame response has the undesirable feature of depending on the flame location 
inside the computational domain. Nevertheless, we feel that the present simulations 
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120 125 130 135 140 145 

Time/ Flame time 

FIGURE 6. Case 1 with ( Z ) = 0.4, A Z = 0.2 (at inlet), and Iz/^i — 11, where 

S° L is the thermal thickness of the perfectly premixed flame with Z = {Z). : 

reduced flame thickness ; ; reduced maximum temperature O' ; + : reduced 

flame speed r* ; * : flame distance to inlet. 


can still be used to describe the basic features of partially premixed flames. 

The example presented in Fig. 5 corresponds to a perturbation in Z with alter- 
native fuel rich (Z > Z 9t ) and fuel lean (Z < Z 9t ) pockets. The excess fuel and 
excess oxidizer that are not consumed by the premixed flame will burn in a diffu- 
sion flame. The intensity of this post-diffusion flame is rather low and in the case 
of Fig. 5, some unburnt fuel is found at the outlet of the computational domain. In 
a similar simulation, but with (Z) = 0.45, there is no leakage of fuel. 

Figure 6 presents typical time variations of the different diagnostics used to char- 
acterize the flame response. When the flame meets a pocket with a mixture compo- 
sition close to stoichiometry, the flame speed and temperature increase, the flame 
thickness decreases, and the flame moves upstream in the computational domain. 
The converse is true when a pocket of mixture composition further from stoichiom- 
etry reaches the flame. Variations in flame speeds are large, with 0.5 < r f < 1.5, 
and show deviations from a sinusoidal evolution: the time required for the flame to 
cross a given constant- Z pocket increases as Z moves away from stoichiometric con- 
ditions. This bias accounts for a reduced overall mean combustion rate compared 
to the perfectly premixed case, thus R f < 1. 

Depending on the values of (Z), A Z, and the effect of partial premixing 
on the mean reaction rate can either be positive, with R f > 1, or negative, with 
R! < 1. Figures 7 and 8 show that this effect remains weak, however, except for 
conditions close to the flammability limit. In Fig. 7, mixtures with AZ = 0.2, 
and (Z) below 0.38 are quenched, while they would burn if perfectly premixed 
(AZ = 0.0). Similarly, in Fig. 8 mixtures with AZ = 0.2, (Z) — 0.4, and lz/^\ > 14 
are quenched while they would burn if perfectly premixed. Fig. 8 also shows a 
comparison between the adiabatic and nonadiabatic simulations. Differences are 
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FIGURE 7 . Case 1 with variable (Z). □ : R f (reduced overall mean combustion 
rate) in partially premixed flames (A Z = 0.2 at inlet) ; ■ : R' in perfectly premixed 
flames (AZ = 0.). 



Figure 8. Case 1 with variable /^, (Z) = 0.4, and AZ = 0.2 (at inlet). □ : R 
(reduced overall mean combustion rate) in adiabatic flames ; * : R in non-adiabatic 
flames ; * : asymptotic value S a for adiabatic flames ; o : asymptotic value S a for 
non-adiabatic flames. 

small until transition to extinction is observed in the non-adiabatic case. In both 
cases, as lz becomes very large, the mean flame speed tends to an asymptotic value 
S a given by the following expression: 

S a - — - 2 - - r (19) 

where S^T = Si((Z) — AZ/2) and 5^ = Sl{(Z) + AZ/2). In Fig. 8, without heat 
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loss, S a /S L {Z st ) = 0.51 ; with heat loss, S a = 0. 

In summary, partial premixing in the one-dimensional case leads to strong tempo- 
ral variations of the laminar flame structure, and in particular to strong fluctuations 
in the instantaneous values of the flame speed Si and the mass burning rate m. 
In the absence of quenching, these variations tend to cancel in the mean, and (Si) 
and (m) remain close to the values of Sl and m obtained in perfectly premixed sys- 
tems. However, quenching induced by partial premixing has been observed in the 
case of strong variations in mixture composition, characterized by large amplitudes 
of A Z > 0.2, or large length scales of lzl&\ > 10. 

5. Case 2: Two-dimensional, y-inhomogeneous, steady, laminar flames 

In this configuration, inhomogeneities in the reactant species exist in the direction 
tangent to the flame, allowing the solution to converge to a steady state. This 
configuration is depicted in Figure 2, where the mixture fraction at the inlet is 
given by: 

. . A Z 

Z = (Z)- — cos 

For two dimensional flows, the parameter space becomes larger than the one- 
dimensional flows discussed previously, and we restrict ourselves to varying (Z) 
and A Z while maintaining lz constant. As these parameters are varied, we expect 
both the flame structure and propagation speed to change. As a result, the flame 
can advance or recede out of the computational domain. To avoid this problem, 
the inlet velocity, which remains uniform, is adjusted to accommodate changes in 
the flame speed. This procedure has been used in partially premixed combustion 
(Ruetsch et ai 1995 and Ruetsch and Broadwell 1995) and results in a steady-state 
configuration. This allows a well defined flame speed to be assessed in each run. 
Note that by defining the flame speed as the inlet velocity required to reach a steady 
state, we are considering a displacement speed. 

As in the one-dimensional case, the mixture fraction is greatly modified from 
the time it is specified at the inlet to the time it reaches the flame. The range in 
mixture fraction at the flame surface is affected by several phenomena, including 
diffusion and the strain induced by the flame. Strain does not directly affect the 
mixture fraction, but does so implicitly by modifying the mixture fraction gradient 
in the lateral direction, which alters mass diffusion. The range of mixture fraction 
on the flame surface for all cases is shown in Fig. 9 as a function of the average 
mixture fraction. In addition to the reduction in mixture fraction range, the mini- 
mum and maximum values are no longer centered around the average value of the 
mixture fraction. The reason for this asymmetry becomes clear when we examine 
the structure of the flames when exposed to gradients in the mixture fraction. 

5.1 Flame structure 

The reaction rates and streamlines for flames subjected to different levels of (Z) 
are displayed in Fig. 10. For (Z) = Zst = 0.5, we observe two leading edge flames 
within the domain. Since (Z) is at the stoichiometric value, we expect two equidis- 
tant leading edge flames and two equidistant troughs. As we decrease (Z) from the 
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(Z) 

FIGURE 9. Range in mixture fraction on flame surface as a function of (Z). The 
different symbols correspond to different values of AZ at the inlet according to 
the following: • represents the homogeneous case (AZ = 0 at inlet), + represents 
AZ — 0.2, o represents AZ = 0.4, and A represents AZ = 0.8. For the chemical 
scheme used in this case, Zst = 0.5. 

stoichiometric value, this symmetry no longer exists. For the case with (Z) — 0.45, 
we still have two stoichiometric points on the flame surface, although they have 
moved closer together. For the other cases of (Z) = 0.4 and 0.35, stoichiometric 
points no longer exist on the flame surface. In these cases, the leading edge is 
located where the mixture fraction is closest to the stoichiometric value. 

The reason for the asymmetric nature of the minimum and maximum values of Z 
on the flame surface, as observed in Fig. 9, can be easily understood from the flame 
shapes in Fig. 10. Diffusion of species has a longer time to act before reaching the 
flame surface the farther the flame is from the inlet. Therefore, the difference in 
mixture fraction along a horizontal line between the flame’s leading edge and inlet 
is smaller than this difference along a line passing through the flame trough. For 
the case of (Z) = 0.5, the maximum and minimum values of Z are both located 
in the troughs which occur at the same horizontal location, and we have symmetry 
in minimum and maximum values. As we depart from average stoichiometry, with 
(Z) < Zst, the trough with rich composition moves forward and the lean trough 
backwards, so that diffusion has less time to act in the rich branch as compared to 
the lean branch. Therefore, the mixture fraction in the lean branch moves closer to 
stoichiometry. 

Another factor that affects Z on the flame surface concerns the role strain plays 
on species diffusion. The divergence of streamlines in front of the leading edge 
reduces the mixture fraction gradient along the flame surface at that location, thus 
inhibiting diffusion. The opposite occurs in the flame trough, where the gradient in 
mixture fraction steepens due to the convergence of streamlines, accentuating the 
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FIGURE 10. Contour plots or streamlines and reaction rates for simulations with 
A Z — 0.4 and: ( Z ) = 0.5 top left, ( Z ) = 0.45 top right. (Z) = 0.4 bottom left, and 
(Z) — 0,35 bottom right. 



Z Z 


FIGURE 11. Propagation speed as a function of mixture fraction. The speeds are 
normalized by the homogeneous case at stoichiometric conditions, S^Zst)^ on the 
left, and by the homogeneous case at the average mixture fraction, S£({Z) ) on the 
right. In addition to displaying the average mixture fraction of the run with the 
symbols, the range of mixture fraction on the flame surface is shown by the lines 
through each symbol. The legend for the symbols is provided in Fig. 9. 
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FIGURE 12. Flame wrinkling for the simulations. See Fig. 9 for a description of 
the symbols. 


diffusion process. 

As we progress towards lean mixture fractions to the point where stoichiomet- 
ric conditions do not exist on the flame surface, the flame shape changes to the 
point where the spatial extent doubles, as the trough corresponding to rich mixture 
fractions disappears. This, in effect, alters the parameter lz without changing the 
computational domain. This doubling in lateral dimension has an effect on the 
flame speed as well as the flame shape, as is discussed in the next section. 

5.2 Flame speed 

When discussing the change in flame speed due to the inhomogeneous medium, it 
is useful to relate this displacement speed to that of the homogeneous case at both 
the average and stoichiometric mixture fractions, as shown in Fig. 11. In these 
figures, both the average mixture fraction at the inlet and the range of mixture 
fraction on the flame surface are shown by the symbols and lines, respectively. 
We begin discussion of the flame speed examining what occurs when the average 
composition is stoichiometric. Independent of the range in mixture fraction at the 
flame surface, the propagation speed remains that of the homogeneous case. This 
behavior was previously observed (Ruetsch and Broadwell 1995) when studying 
confined flames. For this value of the lateral divergence of streamlines due to 
heat release is greatly inhibited by the confinement, and therefore the heat release 
mechanism responsible for enhanced flame speeds, as in the case of triple flames 
(Ruetsch and Broadwell 1995), is absent. As we depart from stoichiometry in the 
mean, the flame shape changes, effectively doubling /^, and we quickly move into a 
regime where streamline divergence is much stronger in front of the leading edge. 
An increase in flame speed is observed relative to Sl((Z)), and in some cases even 
relative to S° l (Zst)- It is interesting to note that flames with lean compositions 
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FIGURE 13. Average local reaction rate along flame surface normalized by the 
homogeneous case at Zst (left) and (Z) (right), as a function of the average mixture 
fraction. 


along the entire length of the flame, designated in Fig. 11 by lines that do not 
cross Z — Zst = 0.5, can achieve flame speeds greater than the homogeneous 
stoichiometric case. 

As (Z) further decreases, the reduction in reaction rate along the flame intensifies, 
and in spite of the streamline divergence the flame speed drops. As the flame 
speed drops, and along with it the inlet velocity in order to stabilize the flame in 
the computational domain, the mixture has a longer time to laterally diffuse as it 
approaches the flame. It is for this reason that the small range in mixture fraction 
at the flame surface is observed for very lean mixtures, apparent from Fig. 9. and 
is why the flame speed collapses to the homogeneous case. 

Thus far we have concentrated on variations of flame speed with the average 
mixture fraction. We now turn our attention to how the fluctuation in mixture 
fraction about the mean affects the flame speed. As we have already mentioned, at 
mean stoichiometry the degree of inhomogeneitv plays no role in flame speed. As 
the average mixture becomes lean, the flame speed increases as long as the com- 
position at the leading edge stays near the stoichiometric value. For flames where 
the composition along the surface is always lean, the greatest speeds in absolute 
terms occur when the range in mixture fraction is the largest. This feature can 
be explained if we re-examine the flame structure. The streamline divergence de- 
pends on the flame curvature, which itself is determined by the local burning rate 
hence species composition. Therefore, the greatest range in mixture fraction along 
the surface would generate the greatest streamline divergence and increase in flame 
speed. This does not hold when the composition along the flame surface crosses 
stoichiometric values. 


5.3 Fuel consumption 

Having discussed the flame structure and propagation, we now turn our attention 
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to fuel consumption. Due to mass conservation, the global consumption rates given 
by R and R! are equivalent to the flame speed ratios shown in Fig. 11. There are 
slight discrepancies between the global reaction rate and flame speed ratios resulting 
from excess fuel leaving the domain in cases which have stoichiometric values on the 
flame surface. However, the domain is large enough, with a grid of N x = 361 and 
N y = 121, that almost all of the fuel is burned in either the premixed or diffusion 
modes before the flow exists the domain. We therefore use the flame speed ratios 
in Fig. 11 as R and R ' in the following discussion. 

For laminar flames we decompose the global burning rate in terms of the flame 
wrinkling, W, and the average burning rate along the flame surface, { r)s or (r')s, 
according to the following relations: 

R = <r) 5 W; R! = (r')sW. (20) 

The frame wrinkling is given in Fig. 12, and the mean reaction rates along the 
flame surface in Fig. 13. It is clear from Figs. 12 and 13 that flame wrinkling is the 
predominant factor in the global reaction rate modification. 

The predominance of flame wrinkling over reaction zone modification is apparent 
for this steady-state configuration. We must now turn our attention to assessing 
whether this trend prevails when we consider flows with unsteadiness in both the 
scalar and flow fields. We address this issue for unsteady scalar fields in the following 
case, followed by a fully turbulent configuration. 

6. Case 3: 2D, xy-inhomogeneous, unsteady, laminar flames 

Case 3 includes two slightly different configurations, shown in Fig. 14. In Case 
3a, the perturbations in Z correspond to an isolated pair of fuel lean and fuel rich 
pockets, whereas in Case 3b, the perturbations in Z correspond to an infinite array 
of such pockets. Case 3a provides basic information on the impulse response of a 
laminar flame subjected both to normal and tangential Z gradients, while Case 3b 
provides information on the response of a flame to periodic Z forcing. Table II gives 
the run parameters for the different simulations. 

Table II. Simulation parameters for Case 3 


Run 

Case 

(Z) 

AZ 

h/6° L ((Z)) 

h 

N r X Ny 

A 

3a 

0.4 

0.2 

5.0 

0 . 

127 x 127 

B 

3a 

0.35 

0.2 

3.5 

0 . 

127 x 127 

C 

3b 

0.4 

0.2 

5. 

0 . 

127 x 127 

D 

3b 

0.35 

0.2 

3.5 

0 . 

127 x 127 

E 

3b 

0.3 

0.2 

1.9 

0 . 

127 x 127 

F 

3a 

0.4 

0.2 

5.0 

0.031 

127 x 127 

G 

3a 

0.4 

0.2 

6.2 

0.031 

127 x 127 
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FIGURE 14. Case 3: Two-dimensional, ary- inhomogeneous, unsteady, laminar 
flames, (a) an isolated pair of lean ( Z — { Z ) — AZ/2) and rich ( Z = (Z) + AZ/2) 
pockets, of size lz ; (b) an infinite array of such pockets. 


Figure 15 presents a typical snapshot of isocontours of cj/r, Z, c, and Yp , as 
obtained in run A. Note that the generalized reaction progress variable defined in 
(7) is a good marker of the premixed flame front. Because it is affected by mixing 
within the burnt gases, Yp is not a good choice to track the flame front. At the top 
of the figure, the flame is seen to interact with a fuel lean pocket (Z = 0.3), during 
which it decelerates and is convected downstream. At the bottom of the figure, 
the flame crosses a fuel rich pocket (Z = Z 3 t = 0.5), accelerates, and is convected 
upstream. These variations in the local flame displacement speed w correspond to 
strong variations in the local flame structure, as observed in Case 1. They also 
correspond to flame surface production. 

As done in the steady-state situation of Case 2, we now compare the relative 
weight of the two basic effects of partial premixing as indicated in Eq. (20): the 
modification of the flame structure through (r')s, and the generation of flame sur- 
face due to wrinkling, W . Figure 16 compares the temporal evolution of the relative 
contributions of these two terms to the global reaction rate from data obtained in 
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FIGURE 15. Case 3, run A, with (Z) = 0.4, A Z = 0.2, and lzi&\ — 5- Top 

figure, isocontours of: reaction rate dip ( ) and mixture fraction Z ( ). 

Bottom figure, isocontours of: reaction progress variable c ( ) and product 

mass fraction Yp ( ). 



FIGURE 16. Case 3a, run A. Time evolution of the reduced global reaction rate R ' 

( ), the reduced flame surface area W (• ), and the reduced surface- averaged 

mass burning rate (r f )s ( D )• Time is made non-dimensional by the laminar flame 
time 6° l ((Z))/S l ((Z)). 
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FIGURE 17. Case 3b, run C. Time evolution of the reduced global reaction rate R! 

( ), the reduced flame surface area W (• ), and the reduced surface- averaged 

mass burning rate (r f )s ( D )■ Time is made nondimensional by the laminar flame 
time 6 0 l ({Z))/Sl((Z)). 

0.2 


0.1 


0.0 
- 0.1 
- 0.2 
- 0.3 

FIGURE 18. Case 3a, run A. Time evolution of the surface-averaged flame stretch 
(k) s ( ), and its two components: the surface- averaged strain rate (cit)s 

( ) and the surface-averaged propagation term (u>V.n) 5 ( ). Time is 

made nondimensional by the laminar flame time 6° L ((Z) )/ S l{(Z)). 

run A. Data from run A indicate behavior that is similar to the steady-state situ- 
ation in Case 2. Partial premixing increases the global reaction rate, /?'>!, and 
the magnitude of the increase is typically 30-40%. The dominant effect of partial 
premixing is a production of flame surface area, ( r')s ~ 1 and R' % W. 

Figure 17 presents similar results for run C. After an initial transient phase, the 
flame response reaches a limit cycle with periodic time- variations. At the limit 
cycle, the global reaction rate is increased compared to the perfectly premixed 
configuration, R f > 1 ; the magnitude of that increase is small, typically 10%; and 
this increase is related to flame surface production resulting from partial premixing, 
R 9 % W. 

As indicated by Eq. (13), the production of flame surface area is measured by 
flame stretch, and the surface- averaged flame stretch (k) s can be decomposed into 
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FIGURE 19. Test of Eq. (21): I\ pp vs (Aw/Iz)(S°l/ Sl((Z))). o Case 3a (runs A - 
B) ; • Case 3b (runs C - E); ■ Case 3a with heat losses (runs F - G); o A case with 
a single lean pocket. 


a strain rate term, (ar)s = (V.tt — nn : Vu)s, and a propagation term, (u>V.n) 5 . 
Figure 18 presents the temporal evolution of these two components of flame stretch 
and shows that partially premixed effects on stretch are not limited to the propaga- 
tion term. A strong positive contribution of (ar )5 is also observed. This contribu- 
tion corresponds to a modification of the flow streamlines upstream of the curved 
flame, as observed in Fig. 10 for the steady configuration of Case 2. 

It remains, however, that while the details of the temporal variations of ( k) s 
depend on the effects of both hydrodynamic straining and flame propagation, the 
basic driving mechanism for flame surface production is the variation of the flame 
propagation speed w with mixture composition. A simple estimate of the global 
flame stretch induced by partial premixing may then be expressed as follows: 


Kpp 


Ate 

lz 


( 21 ) 


where Ate is the amplitude of the variations of te measured at the flame location 
(due to molecular diffusion and unsteady effects, Ate is somewhat smaller than 
Sl((Z) + AZ/2) — Sl{(Z) — AZ/2)), and lz is the size of the pocket. We use 
the peak value of (k) s observed in the flame’s response to perturbations in Z to 
estimate the global flame stretch. In nondimensional form, we get the following 
estimate for a Karlovitz number induced by partial premixing: 

PP Aw SjXjZ ) ) 

a lz S L {(Z)) 1 ’ 


This relation is tested in Fig. 19 and is found to be satisfactory. Note, however, 
that the values of this Karlovitz number remain small, I\ pp < 0.2. 

In summary, partial premixing in Case 3 leads to both modification of the flame 
structure and production of flame surface area. In the absence of quenching, the 
dominant effect on the mean reaction rate is flame surface wrinkling, R* % W. It 
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Mixture fraction Z 

FIGURE 20. Case 4A. Joint probability density function of the reduced mass 
burning rate, r* , and the flame mixture fraction, Z. Time — 4 l t /u\ 

is always positive, R! > 1, but the magnitude of that effect as measured by an 
estimate of the partially premixed Karlovitz number remains small, K^ p <0.2. 

Table III. Initial conditions for Case 4 simulations 


Case 

( 4 >) 

(Z) 


Z' 

iz/$L(Z at ) 

u'/S L 

lt/S° L (Z 3t ) 

i?e x 

4A 

0.8 

0.049 

0.3 

0.018 

2 

7.5 

2 

75 

4B 

0.8 

0.049 

0.3 

0.018 

2 

2.5 

2 

25 


7. Case 4: 3D, xyr-inhomogeneous, unsteady, turbulent flame 

In this case, partial premixing effects are compared to those due to the turbulent 
motions. The numerical configuration corresponds to a premixed flame propagating 
into three-dimensional, decaying, isotropic turbulent flow, with variable equivalence 
ratio. We refer the reader to Trouve &; Poinsot (1994) for more information on the 
configuration, as well as the initial and boundary conditions. The new feature in 
the present simulations lies in the initialization of the scalar field in the flow of fresh 
reactants: 1/r, Yo, and Ypj 2 are specified according to a model energy spectrum, 
as proposed by Eswaran & Pope (1988). The initial probability distribution of 
equivalence ratio is a pdf with two peaks at <f> = 0.5 and <$> = 1.1. Because of 
turbulent mixing, this distribution quickly evolves to a Guassian pdf. Two different 
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FIGURE 21. Case 4A. Time evolution of the reduced total reaction rate, R f 

( ), the reduced mean mass burning rate, r' (o ), and the reduced total flame 

surface area, ( W ) ( ). Time is made non-dimensional by the initial, turbulent 

eddy turnover time, l t /u f . 



FIGURE 22. Case 4B. Time evolution of the reduced total reaction rate, R f 
(■ ), the reduced mean mass burning rate, r' (o ), and the reduced total flame 

surface area, (W) ( ). Time is made non-dimensional by the initial, turbulent 

eddy turnover time, lt/u r . 

simulations were performed. The run parameters are given in Table III. In this 
table lz designates the integral length scale of the scalar field, u 9 the turbulent 
rms velocity, l t the integral length scale of the velocity field, and Re t the turbulent 
Reynolds number (based on u ' and /*). Cases 4 A and 4B correspond to strongly 
and moderately turbulent flames, respectively. Also, the present simulations use 
the single step reaction mechanism proposed by Westbrook & Dryer (1981). 

In the simulations, partial premixing results in strong spatial variations of the 
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FIGURE 23. Case 4A. Probability density function of flame stretch, k. Stretch is 
made non-dimensional by the laminar flame time S° L ((Z))/ Sl((Z)). Time = 4 l t /u f . 

local combustion intensity along the turbulent flame front, consistent with the find- 
ings from the previous cases. In Fig. 20, this intensity is quantified by the reduced 
mass burning rate per unit flame surface area, r', where r' is seen to vary between 
0.5 and 1.5. r* is also seen to correlate strongly with the local mixture composition, 
as measured by the flame mixture fraction. Interestingly, the correlation is approx- 
imately linear, so that departures of the mass burning rate m from the reference 
value m((Z)) (obtained from a homogeneous, planar, laminar flame) tend to cancel 
in the mean when averaged over the whole flame. This tendency is confirmed in 
Figs. 21 and 22, which present the temporal evolution of the two components of the 
total reaction rate, written for the turbulent case as: 


R ' = ?{W). 


In Cases 4A and 4B, the mean mass burning rate remains within 10% of unity, 
so r' % 1, and the total reaction rate is approximately proportional to the flame 
surface area, R! ^ (If). 

There are two mechanisms responsible for the production of flame area in these 
turbulent simulations: the interaction of the turbulent velocity field with the flame 
surface, and the partial premixing mechanism described in Cases 2 and 3. Eq. (21) 
can be used to determine the relative weight of these two mechanisms. The following 
nondimensional number gives an estimate of the ratio of stretch resulting from 
partially premixing to stretch due to the turbulent motion: 

_A A. 6\((Z)\ I, SMZ)) 

T ~ h «' S L ((Z)) h i J((Z)) »' 

where the turbulent stretch is estimated using the integral time scale of the tur- 
bulence. If l t % lz , Nt may be further estimated as (Z ! / {Z)){Silu t ). Hence, Nt 
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scales as the inverse of the ratio of a characteristic turbulent flow velocity divided by 
a laminar flame velocity. Nt is likely to remain small in most practical situations. 
At the initial time, Nt ~ 0.03 in Case 4A; and Nt ~ 0.1 in Case 4B. 

This last point is illustrated in Fig. 23. Figure 23 presents a typical probability 
distribution for flame stretch, as obtained in Case 4A. Stretch is normalized in 
Fig. 23 by a laminar flame time so that stretch values can be directly interpreted 
as values of the flame Karlovitz number, I\ a . The simulation values of I\ a range 
from -8 to 4. These values are quite large and the simulated flame is beyond the 
domain of possible stretch resulting from partial premixing, (A' pp < 0.2). Similar 
results are obtained in Case 4B. 

In summary, partial premixing in the turbulent case leads to strong variations in 
the local flame mass burning rate, but these variations tend to average out, r' ss 1. 
Due to the much larger values of turbulent stretch compared to partial premixing 
induced stretch, the production of flame surface area by partial premixing remains 
negligible, Nt <0.1. 

8. Conclusions 

Direct numerical simulations of premixed flames propagating into laminar or tur- 
bulent flow, with variable equivalence ratio, are used in this paper to study the 
effects of partial premixing on the mean reaction rate. The flamelet theory is shown 
to provide a convenient framework to describe partially premixed flames. 

Partial premixing leads to strong variations of the local flamelet structure, and in 
particular to strong variations of the mass burning rate per unit flame surface area, 
in. In the absence of quenching, these variations tend to average out and the effect of 
partial premixing on the mean flamelet structure remains limited, { m)s ~ r mX{Z ) ). 
Note, however, that quenching induced by partial premixing has been observed in 
the present simulations, in the case of strong variations in mixture composition, 
characterized by large amplitudes (AZ > 0.2) or large length scales ( lz > 10). 

Partial premixing induces flame stretch and, in the absence of quenching this 
effect, is dominant for laminar flames. It is always positive and will result, in the 
laminar case, in a partially premixed flame burning faster than the corresponding 
perfectly premixed flame. The magnitude of the effect of partial premixing on flame 
surface production is measured by Eq. (22). Typical values of the flame Karlovitz 
number are below 0.2, and this effect will be negligible in highly turbulent flames. 
This has been observed in the turbulent flames of this study, where wrinkling effects 
from partial premixing are small compared to wrinkling created by the fluid motion 
for the given initial conditions. 
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